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The two-dimensional carbon allotrope graphene has recently attracted a lot of attention from 
researchers in the disciplines of Lattice Field Theory, Lattice QCD and Monte Carlo calculations. 
This interest has been prompted by several remarkable properties of the conduction electrons 
in graphene. For instance, the conical band structure of graphene at low energies is strongly 
reminiscent of relativistic Dirac fermions. Also, due the low Fermi velocity of Vp ~ c/300, 
where c is the speed of light in vacuum, the physics of the conduction electrons in graphene 
is qualitatively similar to Quantum Electrodynamics in a strongly coupled regime. In turn, this 
opens up the prospect of the experimental realization of gapped, strongly correlated states in 
the electronic phase diagram of graphene. Here, we review the experimental and theoretical 
motivations for Lattice Field Theory studies of graphene, and describe the directions that such 
research is likely to progress in during the next few years. We also give a brief overview of the 
two main lattice theories of graphene, the hexagonal Hubbard theory and the low-energy Dirac 
theory. Finally, we describe the prospect of extracting response functions, such as the electric 
conductivity, using Lattice Field Theory calculations. 
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1. Introduction 

Graphene is a novel two-dimensional carbon nanomaterial with unusual electronic properties. 
These deviate from the standard theory of conduction due to the linear dispersion of the charge 
carriers In spite of the "relativistic" energy-momentum relation, the electrons in graphene move 
non-relativistically. This is evidenced by the low Fermi velocity of Vp ~ c/300, where c denotes 
the speed of light in vacuum. While the band structure of graphene is classified as semimetallic, the 
material conducts electricity better than silver and outperforms the carrier mobility in sihcon due 
to its high resistance to impurities and chemical damage. Most importantly, graphene is inspiring 
technological innovations and improvements, such as graphene-based transistors and integrated 
circuits, transparent conducting electrodes, solar cells and ultracapacitors. The discovery of this 
revolutionary material in 2004 was awarded the Nobel Prize in physics in 201 1. 

Compared to Quantum Electrodynamics (QED), the electron-electron coupling in graphene 
is enhanced by a factor of c/vp, and is thus roughly 300 times larger than the fine-structure constant 
a ~ 1/137 1^. Moreover, as the screening length diverges at the neutral point, the Landau picture 
which admits the mapping of a strongly interacting electronic liquid onto a gas of non-interacting 
fermions breaks down. In this situation, the electrical conductivity is of particular interest, as it is 
under intense study in suspended graphene, which can be routinely manufactured and electrically 
contacted. The temperature- and frequency dependence of the conductivity is also a central design 
parameter in nanoelectronic applications, such as transistors and current switches. 

Graphene, as a strongly interacting quantum mechanical many-body system, ranks among 
the most challenging problems in theoretical and computational physics to date. Other prominent 
examples where complex behavior emerges in many-body systems at strong coupling include, to 
name a few, high-J^ superconductivity, superfluidity in ultracold atomic gases, and the formation 
of hadrons in QCD. While significant progress has been made in the theoretical sector using an- 
alytical approaches (such as perturbative, renormalization group and self-consistent treatments), 
non-perturbative methods which take full account of quantum mechanical fluctuations offer the 
best prospect for a complete, a priori solution. A great driving force behind our understanding of 
strongly coupled quantum systems is Lattice QCD, whereby the properties of QCD are computed 
from first principles on a discretized space-time lattice. Due to advances in algorithmics as well as 
in computer power. Lattice QCD is able to predict the properties of hadrons (such as the masses 
of baryonic resonances) to percent-level accuracy. The discovery of graphene has also promoted 
the adoption of Lattice QCD methods, in particular the Hybrid Monte Carlo (HMC) algorithm, in 
condensed matter physics ^ ^. In turn, the HMC algorithm has lead to great progress in 
the closely related area of Hubbard-like models in atomic physics, in particular for the strongly 
interacting Fermi gases in the unitarity limit [^. 

2. Graphene at strong electron-electron coupling 

Graphene is often described both as a strongly interacting many-fermion system, and as an 
"ultimate system of non-interacting electrons". Recent experiments with suspended graphene have 
provided increasing evidence for the "missing" electron-electron interaction. While the sponta- 
neous breakdown of the semimetallic phase is yet to be detected, observed phenomena include 
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the fractional quantum Hall effect ||8|] in an external magnetic field, and a strong upward renor- 
malization of the Fermi velocity at the neutral point [^]. Also, a resistive state in a suspended 
graphene bilayer has recently been reported at zero magnetic field |1C], and hypothesized to be 
of interaction-driven origin. Taken together, these developments suggest that an investigation of 
emergent strong-coupling phenomena in graphene is timely, which should also clarify the question 
whether the electron-electron interaction in graphene is key to our understanding of the electronic 
phase diagram and transport, or whether such effects can safely be neglected. The combination 
of Lattice Field Theory and MC calculations, formulated either in terms of a hexagonal Hubbard 
model based on the tight-binding description of graphene, or in terms of a low-energy theory of 
Dirac fermions, is a promising way to proceed as it captures the physics of strong electron-electron 
interaction without uncontrolled approximations. Below is a (by no means exclusive) list of objec- 
tives that are likely to be central in the application of Lattice Field Theory to graphene in the near 
future: 



• The zero-temperature electronic phase diagram of graphene. What is the critical value 
of the electron-electron coupling where the semimetallic properties of graphene break 
down? Is suspended graphene above or below this critical coupling at zero temperature? Is 
it possible to induce an insulating state by a technologically feasible amount of strain (in- 
creasing the lattice constant of graphene by a few tens of percent) to reduce the probability 
for the electrons to tunnel between neighboring carbon atoms? What are the critical expo- 
nents of the semimetal-insulator transition? How is the Fermi velocity v^- renormalized by 
interactions? To what extent is the critical coupling and the associated exponents affected by 
screening of the electron-electron interaction? 

• Finite temperature and gapped phases in graphene. The effects of finite temperature on 
gapped strong-coupling phases in graphene are especially interesting, as the Mermin- Wagner 
theorem prohibits spontaneous breaking of a continuous symmetry in a two-dimensional 
system. However, there is considerable evidence from MC studies by Hands et al. that 
suspended graphene undergoes a Kosterlitz-Thouless (KT) transition into a gapped phase at 
a finite temperature, which is currently poorly known. How is the zero-temperature phase 
diagram of graphene reflected in the finite-temperature properties? What is the critical tem- 
perature of the putative KT transition? 

• The electrical conductivity of graphene. How does the electron-electron interaction affect 
the conductivity o{(0,T) as a function of frequency (O and temperature T? To what extent is 
the minimal DC conductivity o{(0 = 0,T) determined by the electron-electron interaction, 
and how can the measured (see Ref. [^) temperature dependence of the minimal conductivity 
at low T be explained? How close is suspended graphene (and graphene on a dielectric 
substrate) to a quantum critical point where the semimetallic phase breaks down? 

• The viscosity of the electrons in graphene. Theoretical estimates [|ri|] of the ratio of 
shear viscosity to entropy density rj/s have put forward the possibility that the viscosity 
in graphene is unusually small, which was interpreted as an indicator of complex fluid dy- 
namics, observable on length scales of ~ 1 jLtm. What is the magnitude and temperature 
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dependence of the viscosity in the presence of interactions? How does the appearance of 
spontaneously gapped phases affect the viscosity? Within a MC calculation, can sufficient 
statistics be obtained to find a clear signal for the viscosity? 

Superfluidity in a graphene bilayer. In a bilayer with equal populations of electrons and 
holes on the respective layers, the attractive electron-hole interaction induces pair formation 
of electrons and holes. What is the critical temperature for superfluidity of electron-hole 
pairs in a such a graphene bilayer? Depending on the amount of screening, has been 



suggested to be as high as 300 K, while ~ /iK in a more pessimistic scenario [12]. If the 
system is imbalanced (more electrons than holes or vice versa), will decrease. What is the 
critical imbalance where T^. reaches zero? 

• Ultracold fermionic atoms in a hexagonal optical lattice. What similarities exist between 
the electronic properties of graphene and the analogous observables of ultracold fermionic 
atoms (such as potassium-40) confined to a hexagonal optical lattice? The tunneling ampli- 
tude and the interparticle coupling can be controlled over a wide range on an optical lattice, 
the latter using the Feshbach resonance technique Jl3|]. In real graphene, such tuning can 
in principle be achieved using tension and dielectric materials, albeit in a much more cum- 
bersome and limited way. Significantly, the interparticle coupling on optical lattices can be 
increased beyond Ug ~ 2, providing access to physics that may only occur in an "unphysical" 
parameter range in real graphene. 

3. Hubbard theory of graphene 

Graphene can be described as a tight-binding model with (non-local) electron-electron cou- 
plings. The corresponding (non-relativistic) Hamiltonian is 

where ^ annihilates an electron at lattice site i with spin a € {tii} and «. ^= The tun- 

neling amplitudes between nearest- and next-to-nearest neighbor sites are f ~ 2.5 eV and t' ~ 
0.1 eV, respectively. If only the on-site repulsion ?7qq is accounted for, the phase diagram of 



Eq. ( |3.1[ ) has been computed using the Quantum Monte Carlo (QMC) approach of Ref. Q14|]. In 
such a treatment, graphene exhibits a semimetallic phase (SM) at weak coupling, a gapped spin- 
liquid phase (SL) at intermediate couplings and an antiferromagnetic Mott insulator (AFMI) at 
strong coupling. As Uo^^/t ~ 3.3 in suspended graphene, such systems are tantalizingly close to the 
critical coupling U^Q/t ~ 3.5 for a semimetal-insulator transition into the SL phase. 

At least two mechanisms exist which can change this picture. First, J/qq/j increases under the 
application of strain, due to the decreased overlap of the carbon orbitals. A technologically feasible 
amount of strain, increasing the lattice constant of graphene by a few percent, may then suffice to 
induce the SL state. Second, the non-local couplings in Eq. (Oh should not be neglected at low 



densities, where the electron-electron interaction remains unscreened. Indeed, the couplings ?7jj 



have been computed in density functional theory (DFT) by Ref. [15], which yielded U^^/t ~ 2.0 
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with further sizeable contributions at longer ranges. A particularly interesting question is in what 
way such non-local couplings affect the electronic phase diagram, in particular the critical coupling 
for the SL state. Notably, the Hybrid Monte Carlo (HMC) algorithm which enables global 
lattice updates, has recently been successfully employed for systems with a Hamiltonian similar to 
Eq. (pTTl), in particular the Unitary Fermi Gas (UFG). 



4. Dirac theory of graphene 

Relativistic QFT can be applied to graphene, as the electronic dispersion relation becomes 
linear in the vicinity of the neutral point [^. The resulting linearized low-energy theory, valid in 
the vicinity of the "Dirac points" of graphene, can be studied at strong electron-electron interaction 
using the Lattice Monte Carlo (LMC) framework, and complements the abovementioned Hubbard 
approach. As in Lattice QCD and Lattice Field Theory in general, Euclidean space-time is used in 
order to obtain a positive definite probability measure. The resulting low-energy QFT of graphene 
can then be formulated in terms of the action 

I r r 

^E = ^ d^^dt (5,.Ao)^ - £ / cfxdt V/,D[Ao] Wa. ^[^o] = 7o(^o + ''^o) + Vf T^-^,-, (4. 1) 

where i/a (with = V^^7q) is a four-component fermion field, and denotes the Dirac gamma 
matrices. A graphene monolayer corresponds to N^ = 2, while a bilayer is described hy N^ = A 
in the absence of interlayer couplings. Interlayer tunneling can be accounted for by introducing 
two coupled monolayers, giving rise to the characteristic quadratic dispersion of graphene bilayers. 
The electron-electron interaction is mediated by the instantaneous gauge field A,,, with coupling 
constant = e^ /e. Here e is the effective permittivity of the medium (for example a hemisphere of 
Si02 substrate). The strength of the electron-electron interaction is controlled by the "fine structure 
constant" = /{Anevp) (in units where h = c = \), such that low values of ~ 1/300 and 
£ ~ 1 translate into a large electron-electron coupling. In the presence of a dielectric substrate, £ > 1 



due to the increased permittivity. It should be noted that Eq. ( ft. ID is an example of a "reduced gauge 
theory", as the gauge field Aq propagates in one more spatial dimension than the fermion field ^f, 
which is confined to a (2 + 1) -dimensional surface. 



The LMC treatment of graphene is based on a discretized version of Eq. (4.1) on a lattice of 
space-time points. Unlike the QMC treatment of graphene, this lattice has no correspondence to 
the "physical" hexagonal lattice of carbon atoms. Instead, physical predictions are obtained in the 
continuum limit of vanishing lattice spacing. As the continuum limit can only be approached, a 
question arises to what extent the space-time lattice distorts the continuum physics. For this pur- 



pose, different discretizations of Eq. ( fl.l| ) have been developed, which share the same continuum 
limit but emphasize different aspects of the physics on a lattice. In Lattice QCD, a standard choice 
is staggered (or Kogut-Susskind) fermions [|T7[], which only partially respect the symmetries of the 



continuum theory. The use of overlap fermions (see Ref. [18| and references therein), while signif- 
icantly more computationally expensive, allows for a much closer approach to the continuum limit 
in studies of emergent strong-coupling phenomena, such as spontaneous chiral symmetry breaking. 
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5. Computation of transport properties 

A brief outline of the extraction of the transport properties of graphene using LMC calculations 
is given here, in order to provide an idea of the computational intricacies and challenges. An 
experimentally readily accessible and relevant observable is the minimal DC {(O = 0) conductivity 
a(?i = 0, r), which involves the computation of the Euclidean correlator 

G(t)= /j2^(/(T,x)7(0,0)), 7(T,X) = £ l|ir(T,x)7,l/A(T,x), (5.1) 

k=\,2 

where p{(o) is obtained from the correlator G(t) (which is measured from the MC data) as 

^, , rdco^^, . , . , cosh(a)T-a)/2r) 

^(^^ = 70 2^^^-'^)^(-)' ^("'^)= sinh(./2r) • ^^-^^ 

The minimal conductivity is then given by the zero-frequency limit of the spectral function p(«) 
according to 

= lim (5.3) 
where the temperature T is related to the temporal extent A^t- of the Euclidean space-time lattice by 



\/T = aNt, and a is the temporal lattice spacing. The inversion of Eq. (5.2) is complicated by the 
fact that G(t) is obtained by MC calculation at a discrete set of points T,, where the number of data 
points / is typically (^(10) due to CPU power limitations and algorithm scaling, while p{(o) is in 
principle a continuous function of ft). As the integration range in (O is discretized into N^^ ~ ^(10^') 
points, a straightforward inversion is ill-defined. 



The numerical inversion of expressions similar to Eq. ( |5.2| ) has become the focus of intense 
theoretical and computational efforts. Here, we shall focus on the studies at Swansea University 
of the behavior of the conductivity in the vicinity of the deconfinement transition in QCD at high 



temperatures by Aarts et al. in Ref. []19|]. These efforts have led to the development of a promis- 
ing algorithm, which allows for a numerically stable calculation of p{(o) using a combination of 
Bayesian analysis and the maximum entropy method (MEM). In the Bayesian approach, one con- 
structs the "most probable" spectral function by minimizing a conditional probability P[p\DH], 
where D denotes the available data on G{x) and H some additional "prior knowledge". In the 
MEM approach, this additional knowledge is introduced by the "entropy term" 

n( m\ 

(5.4) 



p{(0)-m{(0)-p{(0)ln^^ 
m[(0) 



where X'^ is the standard likelihood function, and a is a parameter which controls the relative 
weight of the data and the prior knowledge, which is introduced through the "default model" func- 
tion m{(o). Unfortunately, the MEM analysis of current-current correlators in Lattice QCD suffers 
from numerical instabilities and poor convergence, which make a model-independent determina- 
tion of a in the limit 0) — > difficult. In this situation, performing the analysis in terms of the new 
quantities 

K{(0,z) = —K{(0,z), p{o),r) = — p{o),t), a^m(o)) = mQ{b + aco). (5.5) 
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was found in Ref. [ p^ ] to circumvent the numerical instabilities at small ft) to a large extent. In 
the context of Lattice QCD, the zero-frequency conductivity has been shown to be independent 
of b to an accuracy of 10 — 20%. While such accuracy in itself is sufficient to produce novel and 
valuable information, it should also be noted that as far as CPU power is concerned, the situation 
for graphene is more favorable than the situation for QCD, due to the lower dimensionality of the 
graphene problem. This advantage can be translated to greater MC statistics, but more significantly 
to a larger extent A^^ in the temporal lattice dimension, which will improve the accuracy and model- 
independence of the MEM analysis. 
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